INFINITE DIMENSIONAL FAMILIES 
OF LOCALLY NONSOLVABLE 
PARTIAL DIFFERENTIAL OPERATORS 

MICHAEL CHRIST, G. E. KARADZHOV, AND DETLEF MULLER 

Abstract. Local solvability is analyzed for natural families of partial differential 
operators having double characteristics. In some families the set of all operators 
that are not locally solvable is shown to have both infinite dimension and infinite 
codimension. 



1. Introduction 

A linear partial differential operator L is said to be locally solvable at a point p if 
there exists an open neighborhood V of p such that for any / e C£°(V) there exists 
u G T>'(V) satisfying Lu — f in V. Whereas there exists a rather complete theory for 
operators of principal type with C°° coefficients [1],[11],[7], the theory for operators 
with multiple characteristics is at present limited to various very special classes of 
examples [4], [9], [10], [8], [2]. 

As a test problem one might study operators 

3 3,k 

where {Xj} is a finite collection of C°° vector fields satisfying the bracket hypothesis of 
Hormander, and each coefficient is C°° and real valued. Examples previously an- 
alyzed suggest that for fixed {Xj}, L should be locally solvable for generic coefficients 
djfi- Most work to date has concerned examples of operators possessing high degrees 
of symmetry and depending on finitely many parameters, but recently the following 
examples have been analyzed [2]: L = X 2 + Y 2 + ia(x)[X, Y] in R 3 = {{x,y,t)}, 
where X = d x and Y = d y + xd t , and the coefficient a depends on the coordinate x 
alone but is otherwise arbitrary. It was shown that for real analytic a, L is locally 
solvable at unless a is a constant function taking one of the values {±1, ±3, ±5 . . . }. 
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This suggests that nonsolvability is far rarer, in relative terms, in the class of oper- 
ators with general variable coefficients than it is within the special classes of highly 
symmetric operators discussed by earlier authors. 

The purpose of the present paper is to exhibit a class of examples in which the 
nonsolvable operators form not a discrete set, but instead an infinite-dimensional 
manifold. 1 Thus nonsolvability is quite not so rare an occurrence as the results of [2] 
suggest it to be. 

In R 2 with coordinates (x,t), define X — d x and Y = x m ~ x dt, where it is always 
assumed that to > 2, and consider 



where a G C°° is real valued. For each real constant b define Cb = X 2 + Y 2 + ib[X, Y], 
and define S m to be the set of all b G R such that Cb is not locally solvable at 0. 
Denote by the n-th derivative of a. Our main result is then as follows. 

Theorem 1.1. Suppose that to > 2, a G C°° is real valued and that L is defined as 
in (1.1). 

• Solvability of L at depends only on the Taylor coefficients of a at 0. 

• Each set S m is nonempty and discrete. 

• If a(0) ^ S m then L is locally solvable at 0. 

• For each even m > 2 and each b G S m there exist polynomials Pf. R i— > R such 
that for all coefficients satisfying a(0) = b, L is not locally solvable at if and 
only if for every j > 1, 



The principal and only novel conclusion is the fourth. Consider the set of all 
equivalence classes of germs of smooth functions a at 0, with a equivalent to b if 
they have the same Taylor expansions at 0. Local solvability of an operator L of 
the form (1.1) depends only on the equivalence class of the coefficient a. The fourth 
conclusion is that for even m there is an infinite dimensional family of equivalence 
classes corresponding to nonsolvable operators. 

Curiously, the situation is quite different for odd to. 

Theorem 1.2. Suppose that m > 3 is odd, a G C°° is real valued and that L is 
defined as in (1.1). 

• If to = 3 then L is locally solvable at unless a(0) G S 3 and a^\0) = for 
every j > 1 . 

1 Our examples are nonetheless locally solvable for generic coefficients. 



L = X 2 + Y 2 +ia(x)[X, Y] 






LOCAL SOLVABILITY 



3 



• If m > 5 is odd then for each b G S m there exist polynomials Py. R i— > R and a 
sei E G N of cardinality at most m — 3 suc/i t/iat /or a// coefficients satisfying 
a(0) = b, L is not locally solvable at i/ and on/?/ i//or every 1 < j £ E, 

a^(0) =P,(a«(0),...a (j - 1) (0)). (1.3) 

There is a simple explanation for why there might be many more nonsolvable 
operators here than in the superficially similar three-dimensional situation studied 
in [2]. If a(0) G S m , then local solvability of L is closely related to the asymptotic 
behavior, as e — > 0, of the unique small eigenvalue A(e) of the ordinary differential 
operator B e = —d 2 + y 2 ^" 1 " 1 ) + (m — l)y m ~ 2 a(ey), as will be seen in Section 2; 
A(e) has an asymptotic expansion in powers of e, and nonsolvability is equivalent to 
the vanishing of all coefficients in this expansion. Thus we have a type of inverse 
spectral problem. The situation in [2] is analogous, except that A depends there on 
two parameters while a still depends on only one, so the inverse problem is formally 
overdetermined and hence it is plausible that A vanishes to infinite order at 00 only 
for very few coefficients a. 

These considerations fail to distinguish, however, between odd and even parameters 
m. When m is even there is an extra symmetry in the problem that is absent for 
odd m: B £ is unitarily equivalent to £_ £ , hence X(e) is an even function. Thus 
the inverse problem is highly underdetermined; the situation is quite analogous to 
the classical inverse spectral problem for — d 2 x + q(x) on [—1, 1] [13]. For odd m the 
symmetry is broken, and the inverse problem has a unique solution for m = 3 and is 
underdetermined by at most finitely many parameters for larger odd m. 

One corollary of our theorem is that the nonsolvable operators are not all reducible, 
via obvious symmetries, to examples with constant or very simple coefficients. 

2. Solvability and eigenvalue asymptotics 

It will always be assumed that m > 2. The notation || • ||, with no subscript, 
denotes always the norm in £ 2 (R) with respect to Lebesgue measure. 

Fix 7] G C^°(R) supported in an open set where a is defined, satisfying n = 1 in 
some smaller neighborhood of 0. Redefine b(x) = [a(x) — a(0)]rj(x) and consider 

L = d 2 x + x^-^dl + i(m - l):r m - 2 (a(0) + b(x))d t . (2.1) 

This agrees with our original operator L in a neighborhood of 0, so it suffices to study 
its local solvability instead. 
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Applying the partial Fourier transform with respect to the coordinate t leads to 
the one parameter family of ordinary differential operators 

A T = d 2 x - ar^-V - (m - l)rx m - 2 (a(0) + b{x)), r E R. 

Lemma 2.1. Suppose there exist c, M E R + and a neighborhood V of stic/i t/iat 
/or a// sufficiently large \r\, 

\\f\\ L * < C\r\ M \\Lf\\ L2 (2.2) 
/or a// / G C^V). T/zen L zs locally solvable at 0. 

This is proved in a more general context in the second paragraph of Section 2 of 
[2] with L replaced by its transpose on the right hand side of (2.2). Replacing 
the transpose by the adjoint results in an equivalent inequality, by Taking complex 
conjugates, and in the present case L equals its adjoint. □ 
Assume that r ^ 0. Substituting x = ey, where e = |r| _1 / m , A T is unitarily 
equivalent either to —r 2 l m B^ or to —r 2 / m B~, depending on the sign of r, where 

Bf = -d 2 y + y 2 ^ ±(m- l)y m - 2 (a(0) + b{ey)). (2.3) 

Each Bf is an essentially selfadjoint 2 operator on L 2 (R), with compact resolvent, 
whose spectrum consists of a discrete sequence of eigenvalues, all of multiplicity one, 
tending to +oo. Denote those eigenvalues by Hq(e) < /if (e) < . . . . Then by Lemma 
2.1, in order to prove local solvability of L it suffices to prove that for both choices 
of the ± sign, 

min |^(e)| > ce M as e -> + , (2.4) 

3 

for some c, M E R + . 

Define S+ to be the set of all a E R such that is an eigenvalue of 

C a = -d 2 +y 2 ^ + (m-l)ay m - 2 . 

Define S m = £+ U (-£+) = {a : a E £+ or - a E £+.}. 

Lemma 2.2. is a nonempty, discrete subset o/R. 

Proof. From the identity 

(C f, f) = / r (|^/| 2 + \y m - l f\ 2 ^ dy (2.5) 

it follows that the domain of C a is independent of a, and hence that a i— > C a is an 
analytic family of unbounded operators in the sense of Kato [6] . Since the eigenvalues 

2 See the proof of Lemma 3 of [3]. 
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of C a are all simple, each eigenvalue therefore extends to an entire holomorphic 
function of a G C. Thus either is discrete, or is an eigenvalue of C a for every 
aeC. 

(2.5) implies that ||£ /|| > c \\f\\ f° r some c > 0, for all / in the domain of £ , since 
the right hand side of (2.5) is bounded below by some positive constant times ||/|| 2 
by an elementary argument. Thus is not an eigenvalue of C , so is discrete. 

C a is essentially selfadjoint for each a G R, and £ is strictly positive. In order 
to show that is an eigenvalue of some C a , it suffices by the intermediate value 
theorem and continuous dependence of the eigenvalues on a to show that there exist 
/ G Cq(R) and (3 G R such that (Cpf, f) < 0. But fixing any / that is supported in 
R + and does not vanish identically, this inequality holds for all sufficiently negative 
P. " □ 

Most of the information we require on the operators Bf is either implicit in [2], or 
follows directly from small modifications of that analysis, often with simplifications. 
Therefore we will not repeat the details of the proofs, 3 but will state the results 
needed and refer the reader to the relevant passages in [2]. 

Lemma 2.3. If a(0) S m then 

min|^(e)| > c> as e -> 0. (2.6) 

3 

Therefore if a(0) S m then L is locally solvable at 0. 

The proof of (2.6) is a simplification of the proofs of Lemma 3.2 and 3.3 of [2]. 
Actually L satisfies maximal and subelliptic estimates in this follows from 

general theory [12], or from elementary arguments as in [3]. □ 
We now restrict attention to Bf; all reasoning will apply equally to B~ after minor 
changes of notation. To simplify the notation we then drop the superscript ±, writing 
for all e G R 

B £ = -d 2 y + y* m -V + (m- l)y m - 2 (a(0) + b{ey)). 

Lemma 2.4. If a(0) G S m then there exists 9 > such that for all sufficiently small 
e G R ; B £ has precisely one eigenvalue in [—9,9]. This eigenvalue tends to as 

Alternatively, Lemmas 2.3 and 2.4 may be proved quite simply by observing that for all small 
e > 0, (1 — e)B — Ce < B e < (1 + e)B + Ce for some fixed positive constant C. This simple 
argument does not work in all cases in [2]. Lemmas 2.5 and 2.7 may be avoided by appealing to the 
theory of globally elliptic operators. 
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For a proof see Lemma 3.4 of [2]. □ 
Assume until further notice that a(0) G S m . Denote by X(e) the unique small 
eigenvalue of B e , and fix a real valued eigenfunction i/j of B satisfying ||-?/>o|| = 1, 
associated to the eigenvalue A(0) = 0. 

Lemma 2.5. There exists 5 > such that 

[ (\d y ^ \ 2 + ^)e s ^dy<oc. (2.7) 

J R, 

Moreover there exists C < oo such that for all < r < 5 and all f G L 2 orthogonal 
to il> , 

( {B^ffe^dyKC ( fe^dy. (2.8) 

See the proofs of Lemmas 3.2 and 3.5 of [2]. □ 
By Bq 1 we mean the following operator. B has a one dimensional kernel and 
cokernel, both spanned by ip . Bq 1 is the unique bounded linear operator on L 2 (R) 
satisfying B$ l g = for any multiple g of ip , B B^ 1 g = g for every g _L ip , and 
B Q l g _L 4> Q for every g G L 2 . 

Lemma 2.6. There exist coefficients Aj G R such that for every N, 

N 

A(e) = ^ A^' + 0(e as e -> 0. (2.9) 

3=0 

See the proof of Lemma 3.7 of [2]. □ 
For our purpose it is necessary to recall from [2] the algebraic formalism determin- 
ing the coefficients Aj. Expand 

oo 

as a formal operator valued Taylor series about e = 0, where 
A> = -d 2 + y 2km - x) + (m - l)y m - 2 a(0) = B 

and 

= y m ~ 2 +i ■ (m - l)a ij) (0)/]\ for all j > 1. 

Writing an eigenfunction with eigenvalue X(e) formally as ip(e) = J2j ¥j e ^ with ip = 
■ipo, then equating coefficients of like powers of e in the formal power series equation 
B £ (p(e) ~ \(e)(p(e), gives the following relation. 

n 

BoVn = ~ £(& " Aj^n-j. (2.10) 
J'=l 
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In order that (2.10) should have a solution in L 2 it is necessary and sufficient that the 
right hand side be orthogonal to ipo, which spans the cokernel of Bq. Since \\ipo\\ — 1> 
this orthogonality condition is equivalent to 

n-l 

A„ = (Woio) + - A j)<Pn-j, V>o> (2.11) 

J'=l 

where the summation on the right is empty for n — 1. If (2.11) is satisfied and the 
right hand side of (2.10) belongs to L 2 there exists a unique solution (p n orthogonal 
to the nullspace of B , and we take always this distinguished solution, so that ip n is 
uniquely determined. (2.11) may thus be restated as 

n— 1 

An = {Mo^o) + £<&¥>n-;, ih) (2.12) 

i=i 

since the terms omitted are zero. 

Once ifi ,...ip n and A ,...A„ have been determined, A n+ i is defined by (2.11). 
Then (2.10) is solvable and determines y? n +i- Thus all (pj and Aj are determined by 
induction on j. Lemma 2.5 and an induction argument demonstrate that for each n, 
(fniu) = 0(exp(— r n |y|)) for some r n > 0, so that the right hand side of (2.10) does 
belong to L 2 (R), and the right hand side of (2.11) is well defined. 

Lemma 2.7. Each function ip n belongs to the Schwartz class S. 

This is proved by induction on n. Thus it is given that the right hand side of (2.10) 
belongs to S, and from Lemma 2.5 we know that ip n and d y ip n decay exponentially, in 
the L 2 sense. Since B has polynomial coefficients, a straightforward bootstrapping 
argument shows that if /, d y f, and B f decay exponentially, in the L 2 sense, then so 
does dyf, for all k. (2.10) and induction on n then imply that each derivative dy(p n 
decays exponentially, in the L 2 sense. □ 

Lemma 2.8. If m is even then for all coefficients a, A n = for every odd n > 1. 

Proof. For even m, B £ is unitarily equivalent to B_ e for all e, as is seen by the change 
of variables y i— > —y. Thus X(e) = X(—e), and hence A n must vanish for all odd 
n. □ 

Lemma 2.9. If a(0) G S m remains fixed then each A n is a polynomial in the quan- 
titities a (1) (0),...a (n) (0). 



The coefficients of this polynomial depend only on m and on a(0). 
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Proof. In order to prove this by induction on n, we claim further that for each 
n there exist finitely many functions (p n j G L 2 {H, exp(r n \y\)dy) and polynomials 
g n ,i[aW(0)...a( n )(0)] such that 

V?n = E ^(« (1) (0), • • • a (n) (0)]w, (2-13) 

i 

for all a with prescribed value a(0). Supposing these claims to hold for all n < N, it 
follows from (2.11) that A^+i has the desired form. 4 By (2.10), 

N+l 

<Pn+i — - X] ^iN+i-j^B^iPj - Aj)(p N+1 _ jti 

j=\ i 

then also takes the desired form. □ 

From (2.10) it is now immediate that for any m,n, if a(0) G S m remains fixed, 
then 

A n = c n • a^(0) - Q„(a«(0) . . . ^(0)) (2.14) 
for certain polynomials Q n which depend also on m and on a(0), with 

c n = [(m - l)/nl] f ^ 2 (y)y m - 2+n dy. 

JR. 

In particular 

c n 7^ for all pairs (to, n) such that m + n is even. 

In the special case m = 2, = {—1, —3, —5, . . . } and all the Q n , c n can in princi- 
ple be computed using Hermite polynomials, since the product of y with any Hermite 
polynomial can be expressed explicitly as a linear combination of two Hermite poly- 
nomials. 

The same analysis applies to the operators B~ , with E+ replaced by — E+. For 
even to, because the function y m ~ 2 is nonnegative, is clearly contained in (— oo, 0), 
hence — E+ n E+ = 0. When a(0) G — S+ we define c n , Q n in terms of _B~ by the 
same procedure as above. When to is odd, E+ = — E+ = E m , since the change 
of variables y h- > — y leaves — <9 2 + y 2(m ~ 2 ) invariant but reverses the sign of y m ~ 2 . 
Thus -B+ is unitarily equivalent to -Bz e , so the corresponding small eigenvalues A ± (e) 
satisfy A + (e) = A~(— e) for all e, and hence the coefficients in their asymptotic 
expansions are related by A+ = (— 1) J A7. We have thus proved the following result. 



Nonlinear dependence stems from the occurrence of ftjtp. 



(2.15) 
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Proposition 2.10. Let m > 2 and suppose that a(0) G S m . Suppose that there 
exists j > 1 s«c/i £na£ m — j is even and 

aV\0)^Q j (a<> 1 \0),...aV- 1 \0)). 

Then 

min (e) | > cleP' as e — > 
/or both choices of the ± sign. In particular, L is locally solvable at 0. 
The case of odd m will be taken up in Section 4. 

3. NONSOLVABILITY 

The purpose of this section is to establish Proposition 3.1, a converse to Lemma 
2.10. For even m we denote by X(e) the small eigenvalue of if a(0) G £+, and of 
if a(0) G -E+ . 

Proposition 3.1. Suppose that a(0) G S m and iaai A(e) = 0(e N ) as e — > 0, /or a// 

JV. Taen L zs not locally solvable at 0. 

We may suppose without loss of generality that a(0) G E+, by replacing a by 
—a if necessary. We are thus given that Aj = for all j > 0, and that formally 
B e (Ee j (pj) = 0(e N ) for all N. 

Set 

3=0 

where A is a large positive integer to be chosen later. Formally -B £ $ £ = 0(e A+1 ). 
Since each <pj is a Schwartz function, and all Aj vanish, we have more precisely 

\d k y B £ <S> £ (y)\<C kM \e\ A (l + \y\r M (3.1) 

for all k, M, as follows from (2.13) and the fact that the coefficient of e j , in the formal 
expansion for B £ & £ about e = 0, vanishes for all j < A. Define 

F T {x) = F^ A \x) = $ r - 1/2 (r 1 / 2 x) 

for t G R + . Given any N < oo, (3.1) asserts that 

\d^A T F T {x)\ < C n t- n {1 + Ir^xl)'" (3.2) 

for all x G R, all r > 1 and all < a < N, provided that A is chosen to be sufficiently 
large. 
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Fix rj G C£°(R) supported in (4" 1 ,4) and = 1 on [2 _1 ,2]. For large A G R + define 
G x (x,t)= I e UT V (r/\)F T (x)dT 

JR. 

, A 
J 3=0 

Since each <pj is a Schwartz function, 

\&ZFt(x)\ < C\ a,2 (l + \ l,2 \x\)~ M 
for any a, M, for all r G (4 _1 A, 4A), uniformly in A. Writing 

d?%G x (x,t) = J e itT (irf V (r/X)d:F T (x) dr, (3.3) 

it follows that 

\d?d2G x (x,t)\ < CA 1+/3+a/2 (l + X 1/2 \x\)~ M (3-4) 

for all a,/3,M, for some C a ^M- Integrating by parts n times with respect to r in 
(3.3) gives 

d?d%G x (x,t) = ±i n+ H~ n J e ltT d^[T^r](r/X)F T {x)] dr. 

Thus 

\d?d°G x (x,t)\ < CX 1+(i+a/2 (l + A 1/2 |x|)- M (A|t|)- ra (3.5) 

for any a, f3, M, m. This bound will be used in the weaker form 

\$d%G x (x,t)\ < C M \~ M whenever \x\ > A" 1/4 or \t\ > A~ 1/2 , (3.6) 

as follows from (3.5) by choosing M,n to be sufficiently large. 

Applying Plancherel's theorem in the coordinate t and recalling that <p = t/j 7^ 0, 
we obtain a lower bound 

IIGaI&cr) > 

for all sufficiently large A, for some c > 0. So for large A there certainly exists 
z = z(X) G R 2 such that |^| < A~ 1/4 and \z 2 \ < A~ 1/2 , satisfying \G x (z)\ > 1. 

By (3.5), the L 2 norm of G x on the region where \x\ > A" 1 / 4 or \t\ > A -1 / 2 is 
0(X~ M ) for all M. By (3.4), VGa = 0(X 2 ) in the supremum norm. Therefore 
\G x {z') - G x {z)\ < 1/2 for all \z' - z\ < A" 3 , for all sufficiently large A. Fix h G 
Cq°(R 2 ) satisfying / h — 1, and set 

h x (w) = h{X\w-z{X)) 
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for w G R 2 . Then h x is supported where \w — z(\)\ < CX 4 A 3 , so 



G x h x > cA- 8 (3.7) 

R2 

for some c > 0, for all sufficiently large A. On the other hand, 

||^a||c* m < C M A 4M 

for any M. 

Fix yet another cutoff function ( G C£°(R 2 ) satisfying ((x,t) = 1 for \(x, t)\ < 2 
and define 

g x (x,t) =G x (x,t)C(\ 1 ^x,\ 1 / 2 t). 
From (3.2) together with (3.6) it follows that 

Lg x (x,t)=0(\' N ) 

in any C B norm as A — > oo, where iV may be made arbitrarily large by choosing A to 
be sufficiently large. Moreover g x is supported in the region where \(x, t)\ < CA -1 / 4 , 
which shrinks to as A — > oo. 

If L were locally solvable at there would exist [5] B < oo such that 



/■ 



gh <B\\h\\ CB \\L*g\\ CB (3.8) 



for any h,g G supported in t)| < -B -1 }- The preceding analysis shows that 
no such inequality is valid. Fix any constant B. Taking g = g x , h = h x and noting 
that L = L*, \\L*g\\ C B = 0(\~ N ) where N may be taken to be arbitrarily large by 
choosing A to be sufficiently large. h x is not also small, but ||/ia||c s ^ CA 4B , so 
the right hand side in (3.8) is 0(\~ N+4B ) = 0(\~ N / 2 ) for any large preassigned N. 
On the other hand, | / g x h x \ > cA -8 . Thus given B, N may be chosen so that (3.8) 
becomes false for all sufficiently large A. Thus L is not locally solvable at 0. □ 

4. The case of odd m 

Suppose that m > 3 is odd, and that a(0) G S m . Then = E~, and we begin 
by discussing B+ and the coefficients Aj = A+. 

With the notation of Lemma 2.6 and of the discussion following it, consider the 
moments 

fj,j = / x^ 2 (x) dx 

for each integer j > 0, where tp = i/j . These depend of course on m and on a(0). 
Define fij = for j < 0. 
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Lemma 4.1. For each j > —2, 

(2m + 2j + 2)/t 2m+j _ 1 + (m + 2j + 2)a/ Wj _ 1 = ±j(j + 1)0' + 2 )^-i- 

(4.1) 

Proof. We have .Bo?/' = 0, that is, 

^" = ( x 2m " 2 + ax m - 2 )*l). (4.2) 

Multiplying both sides of this equation by x^ +1 ip and integrating by parts gives, for 
all j > -1, 

/ (x 2m+ ^ + ax m+ i- v W 2 = -J x 1+1 (^) 2 + j ^'-V 2 , (4.3) 

where all integrals are taken over R with respect to dx. The final integral is not well 
defined when j < 0, but the identity remains valid if this integral is interpreted to 
be = 0. Similarly, for any j > —2, multiplying both sides of (4.2) by x j+2 tp' and 
integrating by parts in all terms gives 

J ((2m + j)x 2m+j - 1 + (m + j)ax m+j ~ l )il) 2 = (j + 2) J x^ 1 ^') 2 , (4.4) 

interpreting the right hand side to be zero when j = —2. Taking an appropriate 
linear combination of (4.3) and (4.4) yields (4.1) for each j > —1. When j = —2, 
(4.4) alone suffices to give (4.1). □ 

Note that fij > for all even j. 

Lemma 4.2. If m = 3 then for every j > 0, fij ^ 0. If m > 3 is odd then there are 
at most m — 3 indices j satisfying /ij — 0. 

Proof. The recurrence relation (4.1) simplifies to 

(2m + 2j + 2)/i 2m+j _ 1 + (m + 2j + 2)a^ m+j _ 1 = for j = -2, -1, 0. 

(4.5) 

^ S m since — d 2 + y 2 *™" 1 ) is a strictly positive operator, so « ^ 0. Since m is odd 
and \ij 7^ for even j > 0, (4.5) then implies that jii ^ for all i G {m — 3,m — 
2, m — 1, 2m — 3, 2m — 2, 2m — 1}. Moreover when m = 3, Hi +m and have opposite 
signs for i — 0, 1, 2 if a > 0, and have the same signs if a < 0. 

First suppose that m = 3 and a < 0. Then /i, > for all < i < 6, and the 
recurrence relation (4.1) then implies that fij > for all j > 0. If m = 3 but a > 0, 
then /ij > for i = 0, 2, 4 and /ij < for i = 1, 3, 5, and the recurrence relation then 
implies by induction on j that fij > for all even j > and jij < for all odd j > 1. 
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For odd m > 3 this reasoning breaks down to a degree, because (4.5) still gives only 
3 initial relations, while (4.1) decomposes into m independent three term recurrence 
relations. The same reasoning as for m = 3 applies to 3 of these, but m — 3 are left. 

Fix some r £ {0, 1,2,... m — 1} and consider the sequence of coefficients 6j = fi r+ i m . 
Then fej > whenever r + im is even, that is, for every second index i. The recurrence 
relation (4.1) is a three term recurrence relation for the sequence bi, for each fixed r. 
We claim that for fixed r, there can be at most one index k for which bk = 0. Indeed, 
if k is such an index, then setting j = km + r + 1 in (4.1) yields 

(2m + 2{km + r + 1) + 2)b k+2 + (m + 2(fcm + r + 1) + 2)a6 fe+ i = 0. 

This analogue of (4.5) may then be used as an initial condition and combined with 
the recurrence relation (4.1) just as in the case m — 3, to deduce that b { ^ for all 
nonnegative i ^ k. □ 

To prove Theorem 1.2 for m — 3, suppose that L is not locally solvable at 0. Then 
a(0) £ S m and all coefficients A+ = (— 1) J A~ vanish. Write Aj = A+. We have 

Ai = (PiMo) = a^(0Vi. 

Since /ii 7^ 0, a^^(0) must vanish. Suppose, by induction on n, that a^(0) has been 
proved to vanish for all 1 < j < n. Thus f3j = for all 1 < j < n, and hence (2.10) 
gives B {p n = 0, so ip n = since it is by definition orthogonal to the nullspace of B . 
(2.11) then simplifies to 

An = (Wo>o> = (m - l)a (n) (0)/i„/n! . 

Since A n is assumed to vanish and /j, n 7^ 0, this forces a^ n ^(0) = 0. 

The reasoning for the case of odd m > 5 is similar, except that if /i n does vanish 
then no constraint can be deduced on a^(0). 
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